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Some properties of d-dimensional disordered models with long-range random hopping amplitudes are in- 
vestigated numerically at criticality. We concentrate on the correlation dimension c?2 (for d — 2) and the 
nearest level spacing distribution P c (s) (for d = 3) in both the weak (b d 2> 1) and the strong (b d <C 1) cou- 
pling regime, where the parameter b~ d plays the role of the coupling constant of the model. It is found that 
(i) the extrapolated values of d-z are of the form c?2 = Cdb d in the strong coupling limit and d-z — d — a c i/b d 
in the case of weak coupling, and (ii) P c (s) has the asymptotic form P c (s) ~ exp( — AdS a ) for s ^> 1, 
with the critical exponent a — 2 — ad/b d for b d S> 1 and a — 1 + Cdb d for b d -C 1. In these cases the 
numerical coefficients Ad, ad and Cd depend only on the dimensionality. 
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1 Introduction 

Quantum phase transitions in disordered electronic systems remain one of the central problems in conden- 
sed-matter physics. Considerable attention has recently focused on both the energy levels and the critical 
eigenfunctions, which strongly fluctuate near the critical point and thus have multifractal scaling properties 
[1-7]. Wave-function statistics can be characterized through the set of generalized fractal dimensions d q 
which are associated with the scaling of the q-th moment of the wave-function intensity. A complete 
knowledge of d q is equivalent to a complete physical characterization of the fractal [1]. 

Among all the fractal dimensions, the correlation dimension c?2 of the eigenfunctions plays the most 
prominent role, since it has been related to the level compressibility [8-11], the spatial dispersion of the 
diffusion coefficient [12] and the anomalous spreading of a wave-packet at the mobility edge [13]. Multi- 
fractality is also related to the probability overlap of eigenstates with energy separation much greater than 
the average separation level. In particular, d,2 describes these density correlations. It has been shown that 
for multifractal eigenstates these correlations decay slowly no matter how sparse those states are [14]. 

Short-range correlations in the energy levels, which are closely related to the localization properties 
of the corresponding wave-functions [15-17], can be described by the critical distribution P c (s) of the 
normalized spacings s. The large s asymptotic behavior of this distribution is still an open question and 
the subject of discussion [18]. 

Metal-insulator transitions (MIT's) depend on the dimensionality and symmetries of the system and can 
occur in both the strong disorder and the weak disorder regime (strong-coupling or weak-coupling regime, 
respectively, of the corresponding field-theoretical description) as well as in the intermediate regime. Each 
regime is characterized by its respective coupling strength [19]. 

The disorder-induced MIT is usually investigated for Hamiltonians with short-range, off-diagonal ma- 
trix elements (e.g., the canonical Anderson model). Other Hamiltonians exhibiting an MIT in arbitrary 
dimension d are those that include long-range hopping terms [20-25]. The effect of long-range hopping 
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on localization was originally considered by Anderson [22] for randomly distributed impurities in d di- 
mensions with the V(r — r') ~ \r — r'|~' 3 hopping interaction. It is known [22,23,25] that all states are 
extended for (3 < d, whereas for (3 > d the states are localized. Thus, the MIT can be studied by varying 
the exponent (3 at fixed disorder strength. At the transition line (3 = d, a real-space renormalization group 
can be constructed for the distribution of couplings [23,24]. These models are the most convenient for 
studying critical properties numerically since the exact quantum critical point is known ([3 — d) and, in 
addition, they allow the low-dimensional cases to be treated, thus using larger system sizes and reducing 
the numerical effort. 

Although a great progress in understanding critical properties of the Id long-range random hopping 
Hamiltonian has recently been made [5, 7, 9, 1 1, 26-37] explicit results for the 2d and 3d systems are still 
lacking. Our aim was to investigate the two previously mentioned important quantities, the correlation 
dimension and nearest level spacing distribution, at the MIT, which have been left unexplored. The last 
quantity could help us to solve the existing controversy about the form of P c {s) at the MIT. 

The paper begins by first giving the model used for the calculations in Sec. |2] The results for the corre- 
lation dimension and for the critical level spacing distribution are presented in Sec. [3]and|4] respectively. 
Finally, Sec.[5]summarizes our findings. 



In order to fully represent the mesoscopic systems we introduce an explicit dependence on dimensionality 
d in the widely studied power-law random banded matrix (PRBM) ensemble [5,7,9, 1 1,26-36] (for closely 
related models see also Ref. [37]). Thus, we consider a generalization to d dimensions of this ensemble. 
The corresponding Hamiltonian, which describes non-interacting electrons on a disordered d-dimensional 
square lattice with random long-range hopping, is represented by real symmetric matrices, whose entries 
are randomly drawn from a normal distribution with zero mean, (Hij ) = 0, and a variance which depends 
on the distance between the lattice sites r ; 



Using field theoretical methods [5,6, 9-1 1,23,25,26], the PRBM model was shown to undergo a sharp 
transition at (3 = d from localized states for (3 > d to delocalized states for (3 < d. This transition shows all 
the key features of the Anderson MIT, such as multifractality of the eigenfunctions and non-trivial spectral 
compressibility at criticality. In what follows, we focus on the critical value (3 = d. The parameter b d in Eq. 
Q is an effective bandwidth that serves as a continuous control parameter over a whole line of criticality, 
i.e, for an exponent equal to d in the hopping elements Hij ~ b d [23]. Furthermore, it determines the 
critical dimensionless conductance in the same way as the dimensionality labels the different Anderson 
transitions. Each regime is characterized by its respective coupling strength, which depends on the ratio 
((|Wji| 2 )/(|Wy I 2 )) 1 / 2 cx b~ d between diagonal disorder and the off-diagonal transition matrix elements 
of the Hamiltonian [19]. 

Many real systems of interest can be described by Hamiltonians Q. Among such systems are optical 
phonons in disordered dielectric materials coupled by electric dipole forces [38], excitations in two-level 
systems in glasses interacting via elastic strain [39], magnetic impurities in metals coupled by an r -3 
Ruderman-Kittel-Kasuya-Yodida interaction [40], and impurity quasiparticle states in two-dimensional 
disordered d-wave superconductors [41]. It also describes a particle moving fast through a lattice of 
Coulomb scatterers with power-law singularity [25], the dynamics of two interacting particles in a Id 
random potential [42] and a quantum chaotic billiard with a non-analytic boundary [43]. 

The two limiting cases of the Id model, b ^> 1 and b <C 1, which correspond to the weak and the 
strong disorder limits, respectively, can be studied via the mapping onto the supermatrix er-model for b 



2 The model 
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large [5, 26] and using the renormalization-group method of Refs. [23, 24] for small b. In particular, one 
finds the following result for the correlation dimension d 2 at the center of the band 



Thus, for the Id version of Hamiltonian Q, which possesses the line of critical points < b < oo, d 2 
changes from to the system dimensionality d = 1 as b increases. Eq. (0 has been numerically confirmed 
by several groups using exact diagonalization [6,27,29,36]. We stress that, unlike the Id PRBM model, 
it has not until now been possible to analytically solve the 2d and 3d disordered models with long-range 
transfer terms. 

The system sizes used are L = 24, 36, 48, 72 and 96 in 2d, and L — 8, 12, 14 and 16 in 3c?, whereas 
b d ranges in the interval 0.01 < b d < 10. We consider a small energy window, containing about 10% 
of the states around the center of the spectral band. The number of random realizations is such that the 
number of critical levels and eigenstates included for each L is roughly 1.3 x 10 6 , except for the larger 
system size in both dimensions, for which this number is about 6 x 10 5 . In order to reduce edge effects, 
periodic boundary conditions are included. The power-law nature of Hamiltonians Q did not allow us to 
use efficient algorithms, such as the Lanzos algorithm, which is usually applied to the study of the MIT in 
the Anderson model, due to the large degree of sparcity of the corresponding Hamiltonian. Instead, we use 
standard diagonalization subroutines. 

3 The correlation dimension 

In this section we investigate critical fluctuations in the eigenfunctions of the two-dimensional model Q 
in terms of multifractal measures, focusing on the correlation dimension d 2 . At the MIT, where the natural 
length scale (the localization length) diverges, strong fluctuations in the wave-functions ?/> M (r) appear on 
all length scales. These fluctuations can be characterized by a set of inverse participation ratios (IPR) [44] 



where d q is a set of generalized fractal dimensions. The index /i labels different eigenfunctions and il 
denotes a 2-dimensional region with linear dimension L. Equation Q is valid for individual states and 
for their ensemble average since the spectrum of multifractal dimensions has universal features for states 
in the vicinity of the MIT [3]. The inverse of 1^(2) roughly equals the number of nonzero wave-function 
components, for which reason it is a widely accepted measure of the extension of the states. Note that in a 
good metal, for which eigenfunctions are ergodic, the IPR scale with size L as I ll (q) oc L~ dl > q ~~ 1 \ whereas 
in an insulator, with localized states, I^(q) oc L°. 

For the computation of r q we used the standard box-counting procedure [3], first dividing the system of 
L 2 sites into Ni = (L/l) 2 boxes of linear size I and determining the box probability of the wave function 
in the i box by Pi(l) = J2r l^v( r )| 2 > where the summation is restricted to sites within that box, and 
ipfiir) denotes the amplitude of an eigenstate with energy e M at site r. The normalized q-th moments of 
this probability constitute a measure. From this, the mass exponents r q (L), which encode generalized 
dimensions d q (L) = T q (L)/(q — 1), can be obtained [45] 




(2) 




(3) 



r,(L) = lim 



i=l 
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Fig. 2 Finite-size corrections to the correlation dimen- 

0.02 



coupling regime (fc> 2 = 0.05) for different system sizes: 
L = 24 (circles), 36 (squares), 48 (diamonds) and 96 
(triangles). The straight lines whose slopes correspond 
to the values of d,2(L) are linear fits to Eq. J4}. 



(squares) and b 2 = 0.05 (circles); solid lines are fits to 
Eq. @. 



where 5 = l/L denotes the ratio of the box sizes and the system size. It should be made clear that the 
calculation of T q (L) is suitable only if the conditions [3] 

a < I < L < £ , (5) 

are satisfied, where £ is the localization or correlation length and a is the lattice spacing (or any microscopic 
length scale of the system). In practice, r q (L) is found by performing a linear regression of In Yli^i Pi (0 
with In S in a finite interval of 6. In order to properly satisfy the previous conditions @ and (0, we take 5 
to be in the interval (0.1, 0.4). Since we are mainly interested in the correlation dimension, we shall restrict 
ourselves to the value q = 2 and so di =t<z. 

Using the exact eigenstates of Hamiltonian Q obtained from numerical diagonalizations, we evaluate, 
for each value of b and L, the numerator on the right-hand side of Eq. Q for decreasing box sizes, and 
then calculate d,2 (L) from the slope of the graph of the numerator vs In 5. Figure 1 provides an example 

of the In S dependence of In p? (I) in the strong-coupling regime (b 2 = 0.05) for different system sizes: 

i=l 

L = 24 (circles), 36 (squares), 48 (diamonds) and 96 (triangles). Clearly, there is no ambiguity in the 
determination of the slopes that correspond to the values of d 2 (L). These slopes are summarized in Fig. 2. 

In Ref. [36] it was shown that the finite-size corrections to T2(L) are of power-law type (~ L^ 1 ) for 
the Id PRBM model. A similar decay was found for the multifractal spectrum f(a) and its singularity 
strength a in d = 1 and 2 [35]. Thus, in order to predict the asymptotic values of d 2 , a curve of the form 

d 2 (L) = d 2 +a 2 /L , (6) 

is proposed. Other forms, such as exponential (~ e _Q2L ) or inverse logarithmic (~ 1/ lna2£) have been 
rejected since we have checked that none of them can adequately describe the d 2 size behavior. 

In Fig. 2, we represent the finite-size corrections for the correlation dimension d 2 (L) of the 2d disor- 
dered system described by Eq. Q in the strong-coupling regime b 2 = 0.02 (squares) and 0.05 (circles). 
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Fig. 3 b 2 dependence of the correlation dimension d-2 
(circles) for the 2d disordered model. Solid lines are fits 
to Eq. corresponding to the limiting cases of weak 
(b 2 3> 1) and strong (fa 2 -C 1) disorder. 



Note that significant finite-size effects are present. The straight lines are linear fits to Eq. © and intercept 
the vertical axis at 0.0894 ± 0.0004 and 0.235 ± 0.002, respectively. 

The extrapolated values of d 2 , as obtained from the previous fits, are shown in Fig. [5] (circles) as a 
function of the disorder parameter b 2 of the PRBM model. They clearly change continuously from as 
b 2 — > to the system dimensionality d = 2 as b 2 — > oo. In the two limiting cases of weak (b 2 ^> 1) and 
strong (b 2 <C 1) disorder regimes d 2 can be well fitted by 



2- 
c 2 b 2 



«2 
b 2 



6 2 »1 
6 2 «1 



(7) 



respectively. These fits are shown as solid lines in Fig. [3] The fitting parameters are a 2 = 0.16 ± 0.06 and 
C2 = 4.76 ± 0.07. Note the similarity of Eqs. @ with the corresponding to the Id case, Eqs. 0. 



4 The level spacing distribution 

This section is devoted to the analysis of short-range energy level correlations in the three-dimensional 
model. The emphasis is on the large s asymptotic behavior of the level spacing distribution at the mobility 
edge. We first review existing analytical results concerning this distribution. 

On the localized side of the transition, states with close energy levels are typically localized at different 
parts of space and have an exponentially small overlap. Their levels are therefore uncorrected and the 
corresponding spacings are distributed according to the Poisson law 



Pp(s) = exp(-s) . 



(8) 



In contrast, in the metallic regime, the large overlap of delocalized states, which are essentially struc- 
tureless, induces correlations in the spectrum, leading to the well known level repulsion effect. If the 
system is invariant under rotation and under time-reversal symmetry (orthogonal symmetry), the normal- 
ized spacings s follow Wigner-Dyson statistics at the infinite system size limit 



Av(s) = -sexp 



(-!•*) ■ 



(9) 



At the disorder-induced MIT, the statistics of energy levels changes drastically and presents distinct fea- 
tures reflecting criticality of the theory. This statistics is characterized by a third universal (i.e., independent 
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of the system size and of the details of the Hamiltonian model) distribution P c (s), which is different from 
both Wigner-Dyson statistics and the Poisson statistics [16,46]. On the one hand, the influence of the 
MIT on the spectral properties was studied in Refs. [16,46] by means of the impurity diagram technique 
combined with scaling assumptions. In these studies, it was conjectured that 



with k « 3.3, the reason for such behavior being that the Thouless energy at the transition point is of the 
order of the average level spacing (e c /A rj 1), and so the levels's repulsion is effective only for s < 1. 

A short-range plasma model with interaction only between closest neighbors [47-50] suggest that 
the universality connected with the spectral fluctuations at the MIT is the intermediate spectral statistics 

P c (s) = 4sexp(-2s). 

On the other hand, by mapping the energy level distribution onto the Gibbs distribution for a classical 
one-dimensional gas with a repulsive pairwise interaction, Ref . [5 1 ] derived the following asymptotic form 



where the coefficient Ad depends only on the dimensionality d, and where the critical exponent a, which 
ranges in the interval 1 < a < 2, is related to the correlation length exponent v and to the dimensionality 
through a = 1 + (dv)^ 1 . 

Neither as regards the numerical description of P c {s) is there any consensus. The exponential decay, 
Eq. JlOfc . of P c (s) has been confirmed by most groups at different MIT's (see Refs. [52-56] and references 
therein), while an exponent a rs 1.2 has been found in Refs. [57,58] from a fit in the whole range of 
spacings to a distribution of the form P c (s) = Bs exp(— As a ) or, indirectly, from the two-point correlation 
function of the density of states [59]. Anyway, the behavior il Q with some nontrivial 1 < a < 2 is what 
one would expect at the mobility edge. 

It should be pointed out that MIT's generically take place at strong disorder (conventional Anderson 
transition, quantum Hall transition, transition in d = 2 for electrons with strong spin-orbit coupling, etc.). 
In this regime, the predicted [5 1 ] exponent a = 1 + (dv) ~ 1 slightly deviates from unity, making it relatively 
difficult to see on the numerically calculated tails of P c (s) (e.g, at the standard Anderson transition in 3d 
a w 1.2). To overcome this problem, it is necessary to investigate transitions which occur at the opposite 
limit (weak coupling regime). This area has recently been investigated for the model Q when d = 1 and 
2 [18], and an exponent a close to the Wigner-Dyson value a = 2 has been found at large values of b and 
b 2 , respectively. 

Here, we make a similar study for the more realistic and interesting case d = 3. From results of detailed 
high precision numerical investigations, we will show unambiguously that Eq. il It is indeed correct while 
the validity of Eq. (II Oi is limited to the case of very strong disorder (strictly at the limit of infinity coupling 
strength). In addition, we find that the exponent a in Eq. il It continuously varies between 1 and 2 as the 
coupling strength of the Hamiltonian model changes from to oo. 

For the computation of P c (s), we unfold the spectrum in each case to a constant density, and rescale it 
so as to have the mean spacing equal to unity. In order to diminish the magnitude of the relative fluctu- 
ations and to analyze the asymptotic behavior in detail, it is more convenient to consider the cumulative 
level spacing distribution function J(s) = f P(s')ds'. Note that the integration does not change the 
asymptotic behavior of P(s). The Wigner surmise, Eq. (|9j, and the Poisson distribution, Eq. (|8}, yield 
Pw(s) — exp(— 7rs 2 /4) and ip(s) = exp(— s), respectively. 

Figure |4]displays our results for the integrated probability I c (s) of the 3d model for L = 14 at 6 3 = 
0.02, 0.05 and 5, which are depicted consecutively from bottom to top. Dotted and dashed lines, which 
correspond to Iw( s ) an d Ip( s )> respectively, are given for comparison. A gradual crossover in the large 
s tail of I c (s) from the Poisson to the Wigner-Dyson limiting forms as one increases the inverse coupling 
constant b 3 of the model can clearly be seen. So, we can therefore expect an exponent a in Eq. (lilt . 



P c (s) ~ exp(-Ks) , s»l, 



(10) 



forP c (s) 



P c {s) ~ ex-p(-A d s a ) , s > 1 , 



(ID 
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Fig. 4 Integrated probability I c (s) of the 3d system for Fig. 5 Log-log plot of the integrated probability I c (s) 
L = 14 at b 3 = 0.02, 0.05 and 5 (from bottom to top), of the 3d system at b 3 = 0.02, 0.05 and 5 (from bottom 
Dotted and dashed lines are Iw(s) and Ip(s), respec- to top) and different system sizes L — 8 (circles), 12 
tively. (squares) and 16 (diamonds). Dotted and dashed lines 

are 7\y(s) and Ip(s), respectively, and the straight lines 

are fits to Eq. <1U . 



which spans the interval [1,2], in agreement with Ref. [51]. For the Id and 2c? cases the behavior is quite 
similar [18]. 

Next we consider the behavior of I c (s) with system size L. The results for s large of the critical I c (s) 
for the 3d model at different values of b 3 are shown in a log-log scale in Fig. [5] for different system 
sizes: L — 8 (circles), 12 (squares) and 16 (diamonds). Note that I c (s) is an L-independent universal 
scale-invariant function that interpolates, as previously mentioned, between Wigner and Poisson limits. 
This result confirms the existence of a critical distribution exactly at the transition. Dotted and dashed 
lines correspond to I\y( s ) an d Ip(s), respectively. We checked that the normalized variances of P c (s) are 
indeed scale-invariant at each critical point studied [60,61]. The straight line behavior of the data in such 
a plot at all values of b 3 considered is undoubtedly consistent with a b 3 dependent exponent a in Eq. (II It . 
The values of b 3 reported are 0.02, 0.05 and 5, from bottom to top. The best fit to Eq. (II Q in the interval 
2-5 < s < 5 for small b 3 and 2.5 < s < 4 for large b 3 , yields a = 1.008, 1.039 and 1.901, respectively, 
thus confirming the result of [51]. Note that for the large energy ranges considered, where I c (s) vary by 
one to three orders of magnitude, the quality of the fits, which are represented as solid straight lines, is 
evident. 

Finally, the disorder dependence of the critical exponent a, as obtained from the previous fits for the 
3d system (circles) is shown in Fig. [6]in the broad range of the parameter b 3 of the 3d PRBM model. It 
clearly changes continuously from the Poisson value a = 1 as b 3 — » to the Wigner-Dyson value a = 2 
as b 3 — » oo. In the two limiting cases of weak (b 3 3> 1) and strong (b 3 <C 1) disorder regimes it can be 
fitted by 

"J* 1 , 02) 
b 3 < 1 

respectively. These fits are shown as solid lines in Fig. [6] The fitting parameters are 0,3 = 0.18 ± 0.04 
and C3 = 0.96 ± 0.07. Note that at b 3 1 the condition |r, — r j\/b S> 1 is not completely fulfilled for 




© 2003 WILEY- VCH Verlag GmbH & Co. KGaA, Weinhcim 



8 



Sh. First Author: Critical properties 



Fig. 6 b s dependence of the critical exponent a (circles) 
for the 3d disordered system. Solid lines are fits to Eq. 
1121 corresponding to the limiting cases of weak (6 3 
1) and strong (b 3 -C 1) disorder. 



"lO" 2 10" 1 , 10° 10 1 

b 

the system sizes considered and for the largest b 3 reported the corresponding a saturate at smaller values 
than predicted by the previous equation. From eq. (I12> . the Poissonian tail of P c (s), eq. ( 1101 . is recovered 
for large spacings at the limit of very strong coupling b 3 — > 0. So, we conclude that in the case of very 
strongly coupled Hamiltonians only, eq. (II It loses its validity and eq. dlOl applies. 

5 Summary 

We have calculated the correlation dimension of the eigenfunctions and the nearest level spacing dis- 
tribution P c (s) of non-interacting electrons on d-dimensional disordered models (d = 2 and 3) with long- 
range random transfer terms at criticality in the whole range of the coupling constant b~ d . The leading 
finite-size corrections to efe decay algebraically with exponents equal to — 1. At the infinite-size limit, it is 
found that d^ is of the form d2 = Cdb d for small b d and d% = d — ad/b d for large b d . P c (s) is found to be 
scale-independent at all values of b~ d . The large s part of P c (s) obtained is shown to have an exp(— AdS a ) 
decay with 1 < a < 2. Finally, we determined the disorder dependence of a in both the strong (b d <C 1) 
and the weak (b d ^> 1) coupling regimes. At the limit of very strong disorder b d — > 0, we found that 
a — ► 1 and so we obtain the expected results of the Poissonian decay predicted in Refs. [16,46]. 

Acknowledgements The author thanks the Spanish DGESIC for financial support through project numbers BFM2003- 
03800 andFIS2004-03117. 
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